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^5^ ' Abstract. In the paper, analysis of a quantum optical system — three-level atom in a 

•*■?»■ ' quantum electromagnetic field is given. Evolution operators are constructed in closed 

^~^ form. 

(N 

^H ■ 1. Introduction 

•^r I In theory of computation it is known that the three-state element is most efficient for 

classical computers, but its realization remained a problem for many years. The progress 

was achieved in 90-ties (see [Hj, 0, |H])- Optimal number of states of the physical system for 

the expression mod p, where p is a prime, is 3. Indeed, if we denote by Np this number, it 

K^ . is known that Np ~ f{p)N2, where f{p) — -^^^ — , we obtain that the function jf- have the 

\mml I minimum when p — 3. The necessary mathematical theory in classical many- valued logic is 

f~^ . given in 9 . 

^^ ' In the quantum world the many-level systems are ordinary ones. Thus in this paper we 

^^ . will consider one possible version of a quantum processor based on a many-level quantum 

system. 

In |l()j the extension of the universal quantum logic to the multi- valued domain is consid- 

1-^ ' ered, where the unit of memory is the qudit. Thereby it is assumed that arbitrary unitary 

SpH. operators on any d-level system, where d > 2 is any number, can be decomposed into logical 

gates that operate on only two levels at a time. For example, the linear ion trap quantum 

computer jTJ uses only two levels in each ion for computing although additional levels can 

r^ ' be accessed, and are typically needed, for processing and reading out the state of the ion 

Q^. (see ^ni and the literature cited there). 

In the case when d = 3, the set of universal gates is built in the following manner. 
Let us introduce the following unitary operators on C^ ^ C^ : 

^ : / 1 \ / Co -C2 

C^ ■ ^3(<^)= 1 ,Z3(C0,C1,C2)= -Ci C2 

V e'"^ / \ -C2 C2 C2 

It can be easily checked that X3{(f)) and ^3(00, ci, C2) act on the states |0 >, |1 >, |2 > as 

X3(0)|O >= |0 >,X3(0)|1 >- |1 >,X3(<^)|2 >= e''^|2 >; 
Z3(0)|O >= co|0 > -C2I2 >,Z3((/))|1 >= -ci|0 > +C2\2 >,Zs.\2 >= -C2IO > -fc2|l > +C2I2 > 
Furthermore, the gate Z^^cq, 01,02) satisfies the following condition: 

f Co -C2 \ 

^3(C0,C1,C2)= -Ci C2 (Co|0>+Ci|l>+C2|2>) = 

\ -C2 C2 C2 / 

= (-C2C0 + coC2,-C2Ci +C1C2, |cop + |ci|^ + |c2p = |2 >). 
In jni it is proved that the set of gates {X3(0), ^3(00, ci, C2), G2[X3((/))], G2[^3(co, c)l, C2)]} 

is universal for a three-level quantum system, where G2 [Y] = r, t^ I and / is the identity 

matrix. 

Below is given another approach of description of quantum system which gives the sin- 
gle qutrit gates, in particular, we obtain in closed form evolution operators depending on 
parameters. 
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2. Atom-Photon Interactions Hamiltonian 

We will build the special case Hamiltonian of a hydrogen-like atom moving and interacting 
with a quantized electromagnetic field, at that using the notations of ^ and [2j- The general 
form of such Hamiltonian, including the center-of-mass motion and all the modes of the 
interacting field, is presented in [I]. For our purposes we consider the following conditions 
on atom-photon interactions: (i) All the atom-photon interactions are the electric dipole; (ii) 
Only three atomic levels are included in the interaction; (iii) Two quantized resonator modes 
are interacting with this three-level system; (iv) Each of these modes interacts with only 
one couple of levels; (v) So, from three possible couples of levels, only two of them interact 
directly. Levels of the third couple interact only by means of an intermediate level; (vi) In 
this work we permit, but do not consider the center-of-mass motion and large detuning from 
resonances, which cause the state entanglement. 

Let us consider atom-photon interactions Hamiltonian in the electric dipole approximation 

(1) ^r.E = -pE(i?, t) = -erE(R, t), 

and transform it for three-level atom interacting with two modes of quantized fields of 
radiation. For these purposes previously we must discuss the following operators. 

Hamiltonian for internal degrees of freedom of non-interacting atom. We con- 
sider internal degrees of freedom of a three-level atom. Due to assumptions (i-vi) only two 
couples of energy states are involved in interaction. In all cases the intermediate state with 
energy Ei = fwji is labeled by the state vector \i >. From other states, the state with 
higher energy is named excited, denoted Ee = fkUe and labeled by the state vector |e >. 
The remaining energy level is named ground level, denoted Eg = hiOg and labeled by the 
state vector \g > (see Fig.l). All of them are eigenstates of the Hamiltonian Hatom of 
non-interacting three-level atom: 

(2) Hatom\e >= fi-^e|e >, Hatom\i >= ^il* >' Hatom\9 >= ^gls > ■ 

Using the completeness relation |e >< e| + |i >< i| + |(? >< g\ = 1, we can write the 
arbitrary atomic operator O as 

(3) 6 = i6i= J2 b><j|oi/></l- 

j,j'=e,g,i 

Utilizing JSJ and the orthonormality condition < j\\j' >— dij' for the Hamiltonian of a 
three-level atom, one can easily obtain 

(4) Hatom = (|e >< e| -|- \i >< i\ + \g >< g|)i?atom(|e >< e| -|- \i >< i\ + \g >< g\) = 

Ee\e >< e\+E^\i >< i\ + Eg\g >< g\, 
where we use ^ and above notations E^^ — fiWei Ei — hiUi, Eg = fkUg. Then 

|e>=j ,<e| = (l,0,0),|z>= 1 , < z| = (0, 1, 0), |g >= | | , < ff| = (0, 0, 1), 

Thus, the Hamiltonian Q) has the following matrix form 

(5) Hatom — -E-e 





The dipole moment operator. Now we can use the completeness relation once more to 
express the position operator f in terms of energy eigenstates. As each state, wave functions 
^j (r) of a free atom have definite parity, then the diagonal matrix elements equal zero 



<i|f|j>=y"|*j(r)prdV = 0. 



Indeed, if the function |^j(r)p is symmetric and the operator r antisymmetric, the expres- 
sion under the integral sign is antisymmetric. For the nondiagonal matrix elements we have 
the following expressions 



(6) e < z|f|5 >= e / **(r)r*g(r)dV = p,g, 

and 



In the same way 



e < e|f|i >= e / **(r)r*g(r)rf 



r3 

^ Peit 



and 

e < i|f|e >= e / *;(r)r*,(r)d3r = p*,, = p,e- 

Thus, the dipole moment operator ef is expressed in the following manner 

(7) ef = pig\i X g\ + pgi\g >< i\ + pet\e >< i\ + pie\i >< e\. 

Note that this operator describes transitions from the ground state \g > to the intermediate 
state |j >, from intermediate to the excited state |e >, and the process in reverse order back 
to \g >. For verification let uss apply the operator ef (9) to the ground state \g > . Then 
we get 

ef|g >= pig\i >< g\\g > +pgi\g >< i\\g > +pe»|e >< i\\g > +pie\i >< e\\g >= pig\i > . 

Similarly 

ef|e >= p,g\i >< g\\e > +pgi\g >< i\\e > +pei\e X i\\e > +pie\i X e||e >= pie\i > 

and finally 

ef|z >= pig\i X g\\i > +pgi\g X i\\i > +pei\e X i\\i > +pie\i >< e\\g >= pg^\g > +pei\e > 

Thus, the operator 

/ \ / 

(8) a,g = \i><gl \t><g\=[ 1 (0,0,1)= 1 

\ / \ 

causes the atom transition from the ground to the intermediate state and can be interpreted 
as a creation operator for the atom in the intermediate state \i >. 
On the contrary, the operator 

/ \ / 

(9) ag, = \gxi\, \gxi\^ [ (0,1,0)= 

V 1/ Vo 1 

annihilates the atom in the intermediate state and so represents the annihilation operator 
for the intermediate state. 
The operator 

/1\ /O 1 

(10) ae^ = \exil \exi\= \ (0,1,0)= 

\ / \ 

causes transition of the atom from the intermediate to the excited state and can be inter- 
preted as a creation operator for the atom in the excited state |e >. 
The operator 

/ \ / 

(11) (Te. = |z >< e|, |i >< e| = 1 (1,0,0)= 1 

\ / \ 

annihilates the atom in the excited state and so represents the annihilation operator for the 
excited state. The dipole moment operator Q is defined in accordance to the expressions 
(|5|)- (|ll|l . of creation and annihilation operators, to be 

(12) ef = pig&ig + pgi&ig + pei^ci + PlC^ie- 
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The electric field operator. Now let us discuss two modes of radiation field. We 
assume for determinacy, that the radiation field is formed in three-mirror standing wave 
resonator and decompose the electric field strength in the resonator as follows 

(13) E(R, t) = ToaMnKa - a+) + .Fo6Uh(R)i(6 - S+), 



where Tqj = ^J^, j = a, b. 

Here Uj (R) represents the mode functions of the resonator at the position R of the atom. 
ria and rift are oscillation frequencies of the field. In a classical treatment, a+, fe+ and a, 
b are time dependent amplitudes. In our present treatment they are operators which obey 
the commutation relations 

(14) aa+ - a+a ^ 1, 66+ - 6+6 = 1. 

With their aid the quantized field Hamiltonian can be written in the form 

(15) ilfieid = Mla'a^a + hribb^b. 
Here we neglect the zero-point energy. 

Back to the interaction Hamiltonian. Now consider the Hamiltonian of electric 
dipole atom-photon interactions again (^ 

HrE = -efE(R,t). 

Substituting in this expression the expressions of the dipole moment operator ef (|12() and 
electric field operator f E(R, t) we obtain 

(16) HrE = -Toai[pigUa{'R)aig + peiUb{K)ffgi]{a - a+)- 

-J^06i[peiUb(R)(Tei + pejUf,(R)CTie] (6 - 6+). 

This expression contains the scalar product pu of the dipole moment and the mode function 
u. As this product is a complex quantity, we represent it in the following manner 

(17) pu=\pu\e'^, 

where ip is the phase. Substituting (|17|l in (|16|l the Hamiltonian of electric dipole atom- 
photon interactions can be written in the following form 

|p»g Ua(R)l ~ i^^ , |peiUb(R)| 

n — "^^^ + — h 



(18) HrE^-n^Oaii ^ CTige'^^'H agie ^"J(a-a^)- 



-;.^o..[^^2^a.e-^ + ^^Sil^^^^e-'-](6 - 6+). 
n h 

Introducing the Rabi frequencies 

t-a^ |p»gUa(R)| ^ ,„. |p«Ufc(R)l ^ 
ga[ti.) — y-Qa, gb[ti.) — T J-Qb 

the interaction Hamiltonian (|i8|l now is given by 

(19) HrE - %a(R)(-j)(a»ge*'^» + a^.e^*'^" ) (fl - a+)+ 

%6(R)(-z)(a«e*'^'' + CT,ee~'^'){b - b+). 
Introducing the difference of phases Sip ^ if a — ^b-, in the Hamiltonian (|i9|l we obtain 

(20) TJrE = %a(R)(-0(<T,<,e''^" + ag,e-*'^°)(a - a+)+ 

%ft(R)(-i)(ae.e^(''^-'^") + a,ee-*(^'^-'^»))(6 - 6+). 
In quantum mechanics only the difference of phases is important, so in H20() we can put 
ipa — 7r/2 and finally obtain 

(21) HrE - %a(R)(a.g - ag,){a - a+) + %ft(R)(<7e.e-^^'^ - a.ee^^'^)(6 - 6+). 

Expression 1)21(1 is the final form of the exact Hamiltonian responsible for all atom-photon 
interactions between the three-level atom and two modes of electromagnetic field. Just this 
Hamiltonian will be used hereinafter. 



Total Hamiltonian. Joining up the above main results Q), IjlSI) . (|21|l . and taking into 
account, that 

CTee = |e >< e|, 0-ji EE |i >< i\, Ogg = \g X g\, 
one can build the total Hamiltonian of the whole atom-field system 

P2 . 

(22) H — —— + H field + Hatom + H^E = 

2M 



2M 



MlaO-^Ci + hfltb^b + Ee(Tee + Ei&u + E„a. 



a'^gg 



+hga{R){&^g - <7gO(a - a+) + ;i.gfc(R)(ae.e-'^'^ - a^eB^'nCb - b+). 

This Hamiltonian includes: the operator of kinetic energy ^rr of the center of mass 
movement; free-field Hamiltonian H field = MlaO^a + Mlhh^b; the Hamiltonian responsible 
for internal states of the atom Hatom = Eeaee + Ei&n -\- Eg&gg and interaction operator 
HrE = %a(R.)('7zg " ^g^)ia " a+) + /^^^(R) (a^^e-*'''^ - &,ee'^'f'){b - b+). Connecting the 
internal states with field operators a — a^ , b—b^ by means of atom- field coupling constants 

5a(R), 5a (R). 



3. Simplification of Atom-Photon Interactions Hamiltonian by Selecting the 
Real Transitions (SRT) and Rotating Wave Approximation (RWA) 

There are two reasons for simplification of the atom-photon interactions Hamiltonian (J2I|I 

(23) Hr-E = figai'R)i^iga^ + ^giO, ~ GigO. - (7iga+) + 

+%b(R)(<7e.&+e*^'^ - a,ebe-'^^ - a.e^e^^'^ - a»e&+e-^^'^). 

Namely, it amounts to confining our analysis to real transitions in which a photon is emitted 
while the atom goes from its upper state to its lower state, or a photon is absorbed while the 
atom goes from its lower state to its upper state. The second reason is mathematical and is 
connected with the averaging procedure - neglecting the fast oscillating terms. Hereinafter 
we will use the interaction representation (Dirac representation) which is defined by the 
internal atomic states and the free-field states and under these assumptions we may reduce 
the Hamiltonian H23|) to the form which is called Rotating Wave Approximation (RWA). 

In a three-level system the simplification procedure of interaction Hamiltonian should 
be carried out separately for each transition scheme. Thus, first of all let us consider these 
transition schemes. In the three-level atom due to (i)-(v) of section 2 there are three different 
configurations possible for transitions between the three levels depending upon the position 
of the intermediate level \i>. Those are shown in Fig. 1 along with the allowed transitions. 
In the configuration of the first diagrams in Fig. 1, the energy of \i > is intermediate between 
the energies of the other two levels. It is called the ladder configuration. Notation 'L' stands 
for this configuration. The energy of |i > in the A or Raman configuration is higher than 
that of the other two levels (second diagram in Fig. 1) whereas the level |i > in the V-system 
is below the other two levels (third diagram in Fig. 1). 



Q, 



a 



Zy 



A. 



K) 



Is) 



a) L- configuration 




a 



k) 




k) 



b) A- configuration 



c) V- configuration 



Figure 1. Three configuration of the three level atom in a two- mode 
field each acting on only one transition. 

3.1. Simplification for L-configuration - SRT. Selection of real transitions. The 

terms agid, aiga~^, dieb and Oeih^ , in the interaction Hamiltonian (|23|l lead to the violation 
of the energy conservation law. Indeed, note that the operators agid and aieb^ annihilate the 
excited atom (in \i > and |e >) and simultaneously annihilate the field excitation (Fig. 2c). 
In the same manner operators digd^ and a^ib'^ create the field excitation while the atom is 
annihilated in \g > {\i >) and created in \i > (|e >) states (Fig. 2d). 
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Figure 2. 

On the contrary, the operators agid^, OiJi^ and digd^ , a^ib either annihilate the excita- 
tion {\i > or |e >) while creating the field excitation (Fig. 2b), or create the excited atom 
(in \i > or |e >) while annihilating the field excitation (Fig. 2a). 

Under the assumptions of real transitions, neglecting the operators o'gid, dieb, digd'^ , 
o'ieb'^ , we may reduce the Hamiltonian (|23ll to the form 



(24) 



TtL n^ TtL 



hgai^R.)i&g^d+ + a^gd) + hg,{Il){&J+e'''^ + aje'^'^- 



3.2. Simplification for A-configuration. Selection of real transitions. For this con- 
figuration the operators cigid, aigd'^ , aieb'^ , cieib of the interaction Hamiltonian lead to the 
violation of the energy conservation law (see Fig.l and Fig. 3c, d). Thus, neglecting these 
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terms in the interaction Hamiltonian (|23|l and keeping UgiO^ ^ criga, Uieh and Ueib^ , and (Fig. 
3a, b), we may reduce the Hamiltonian 123|) to the form 



10:1') 

|g>.k) 



3 

T 



10:1') 

|g>k) 



I'). I') * 

|g).k> -^ 



I'). I') 
|g).k> 



(25) 



Figure 3. 
H^E = Htnt = /J5a(R)(CT<„a+ + a,ga) - %fc(R)(a,e6+e**'^ + ^Je-'^^- 



3.3. Simplification for V-configuration. Selection of real transitions. For this con- 
figuration the operators agia~^, aigd, ai^b, and o'etb^, of the interaction Hamiltonian lead 
to the violation of the energy conservation law (see Fig.l and Fig. 4c, d). Thus, neglecting 
these terms in the interaction Hamiltonian (|23|l and keeping agiO,, digO,, aieb, and &eib, and 
(Fig. 4a, b), we may reduce the Hamiltonian (|23|l to the form 
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Figure 4. 



(26) 



H^E = Hit = -hga{^){ag^a^ 



- cT.ga) + %fc(R)(a,eS+e'*'^ + &e^be-'''P). 



-i5<p\ 



3.4. Interaction representation. Further simplification of the atom-photon interaction 
hamiltonians H21|l - H26|l takes place in the interaction representation, where we can use the 
averaging procedure by fast oscillations. To accomplish this representation we should rep- 
resent the total Hamiltonian of the whole atom-field system (see (I22|) ') 

H — -TTZ + H field + Hatom + ^rE — 



2M 

p2 , , 

— — + Ha + hga(R){aig - (7gi)(a - a+) + hgb(R){d-eie 

by the internal states of the atom and free-field states. Here 

Hq = Hfiid + Hatom = hClaa^a + hQbb^b + Eeffec 

or 

, Wr 
(27) Ho = hflaci+a + Klbb+b + h 



iSip 



-iSif 



(5-6+). 



I- E^ai; 






^g^gg^ 







denotes the sum of free-field H field and free atom Hatom Hamiltonians. 

Now, defining the state vector in the interaction representation |$(^'(i) > via 



(28) 



|$(t) >= exp[-^Hot]\^'^^\t) > 



and substituting H28|) . containing the state vector |<i>(-^)(i) > in the interaction representation, 
in the Schrodinger equation 



We receive 



where 






(29) hI^ = e^p[-Hot]HrE exp[--7?oi] 

n n 

is the interaction Haniiltonian in the interaction picture. Note, that by reason of (vi) of 
section 2, we do not include the kinetic energy operator (corresponding to the center of 
mass movement) in above transformations. Now, let us obtain the exact expression for 
the interaction Hamiltonian H^-J. Substituting Hq (|27|) into the expression (|29l) . under 
conditions, that the field and atom operators commute with each other we obtain 

(30) H^'^ = %a(R)e*^— *(<7,<, - <75,)e~*''°'-*e^''"''^'*(a - a+)e-'f^'"^'''^* 

At first consider the part containing the atom operators. As the matrix corresponding to 
Hatom is diagonal (see (0), HHl), then we obtain the following expression 

(31) exp{-Hatomt) 

n 

Utilizing H3U|) in H31|l with the aid of (|SJ| one can easily obtain that 



,THatomt)aigexp{--^ 

j(uJi-UJg)t 



eXp{-Hatomt)(TigeXp{--Hatomt)a-ig 



= O-igB 

In the same way, taking into account ©, l|10|l and (|ll|l 

eXp{-Hatornt)^gieXp{--Hatomt)o-tg = &gie~''^'^'~'^^'>\ 

expCTHato^t)cj,,e-'''^exp{~lHato„,t) = <7,,e^('^=-'^')*e-^^'^, 

expi^Hatomt)c7,,e'^''exp{-'-Hatomt)a,g = a,;ge-*("=^'^')*e'*'^. 
Thus, the atomic part of the interaction hamiltonian H3U|I is expressed in the following form 

(32) e^^— *(a,g - agOe"^^-""* + e^^— '(ae.e"**'^ - a,ee*^'^)e-^^-<""* = 

Now let us calculate the field operators. The time dependence of the operators in the 
interaction representation is determined by the Heisenberg equations of motion 

(33) — / = -[H field, f] = TiHfiiedf - fH filed), 

where / stands for a, a+, 6 and b. Inserting Hfi^id (I15II in (|33|l we obtain 

(34) jj^\^{f+ff-ff+f). 

Because of the commutation relation 1141) the equation (|34ll is transformed into 

(35) 1/ = -"/"• 
In the same manner 

(36) jj+ = inf+. 



So the time dependance of field operators is defined from H35|) and (|36|l and arc equal 

(37) a{t) = a(0)e-^"°*, ait)+ = a+(0)e*""*, b(t) = b{0)e''^'\ b{t)+ = b+{0)e'^^\ 

Thus, by substituting (|32|l and 1371) into the expression (|30|l we obtain the atom-photon 
interactions Hamiltonian in the interaction representation 

+hga{^R)[ae^e'^'^^-'^''>'e-'^•^ - aiee~'^'^^--'^^'>'e'^'^]{be-'^''' - 6+6^"*). 
Introduce the following notation 

(38) LUi - UJj = LOij . 

Taking into account (|38|l the Hamiltonian is given by the following expression 





hIQ = ;i5,(R)[a,ge'"-* - ag,e-''^-'](ae-'""* - a+e'""* 




+Hgb{^R)[ae^e''^"'e'^^ - a^ee~''^"']{be~'"'* - S+e'""*) 


or 




(39) 




MI) _ 


= %a(R)('Tg,a+e'(""~'^'«)*+a,gae-'("''-'^'«)*-CTg,ae-*(""+^'s)*-e-^ 



+hgbiK){a^eb'^e'^"''~•^"^'+'^'^+ae^be-''■'^''-'^"'^'-'^'^~a,ebe-'^"''+'^^^^^ 

Further transformations require choosing particular transition scheme (Fig. 1). 

L-configuration. For this configuration (Fig.l) 

(40) UJg < LUi < UJe ^ ^ig > 0, Wgj > 0. 

After introducing notations 
and 

Ab = rib — Wei. 

from (|39|) and l|4U|) we obtain the interaction Hamiltonian for L-configuration H^^ 

^rE^^^ = %a(R)(a3.a+e'^°* + a.gae-^^"* - ag^ae"'^*"" +"'■''* - a,3a+e^(""+"-<')*)+ 

-|-%6(R)(6-,e&+e''^'*+''^'^ + (7e»6e-*'^^*-**'^ - ^.J^^-^(^b+^.^)t+^Sv „ ^^^^+g»(n,+a;„)*-*'5v)_ 

It is easy to see that operators (JgiCi^ (JigO^ , (Ji^b and (tJ^^"*" leading to the violation of energy 
conservation law are multiplied by factors containing the sums of field oscillating and atom 
transitions frequencies. Contrary, operators agiO,'^, o'igO,, aieb'^ and aeib are multiplied by 
factors containing the differences (detuning) of aforementioned frequencies A = fi — w. So in 
the interaction Hamiltonian we have two groups of terms. The terms of the first group oscil- 
late approximately with double optical frequencies and their contribution in the Schrodinger 
equation is vanishingly small. The terms of another group oscillate with detuning frequen- 
cies A, and vary slowly. So their contribution in the Schrodinger equation is significant. 
These terms are UgiO^ , cfigO., dieb'^ and Geib. So utilize the averaging procedure - neglecting 
the fast oscillating terms such as Ugicx, UigCL^ , Uipb and Ueib'^ the interaction Hamiltonian 
reduces to the following form 

It is also possible to strongly validate given approach by the mathematical averaging, con- 
trary to heuristic averaging used above. Also the high order contributions can be calculated 
in this case. 



A-configuration. For this configuration (Fig.l) 
After introducing notations 



Uig > 0,U!ei < 0. 



Aa = ila — LOig 

and 

Ab = rib — Wie. 
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we obtain the interaction Haniiltonian for A-configuration -ff^E i-"^ ^^'■^ interaction repre- 
sentation 

Averaging 1)41(1 by fast oscillations leads to 

^-configuration. For this configuration (Fig.l) 

UJg < LU^ < LUe ^ LO.g < 0, LUei > 0. 

After introducing notations 

Aa = ila — OJgi 

and 

Ab = fib — Wei, 

we obtain the interaction Hamiltonian for ^-configuration H!^-^ in the interaction repre- 
sentation 

(42) H^'^''^ = %a(R)(ag,a+e'(^^"+'^-)* + <7,gae-'(""+"-)* - a^.ae-''^"* - a,ga+e'^^' + 

+ngb(R){a^eb+e"^''*+'^'^ + aeibe-'^"^-'^^ - a,Je-'^^''+'^"'>*+'^'^ - ae^b+e-'^^''+'^"^*-'^'^). 
Neglecting in 1421) fast oscillating terms we have 

Therefore, in this section we have shown that in the RWA three-level atom interacting 
with two modes of standing wave resonator, in the interaction representation is described 
by the following Hamiltonians 

L-configuration 
(43) 
^^e"^^ = Hit = hgaiag,a+e-'^'' + a.gae-'^-') + hgbiaj+e'^''+'''^ + a^fee^^^"*"^*'^). 

A-configuration 
(44) 7?(^'^^ = H^^, = hga{iJg,a+e'^-'+a,gae~'^^*) + hgb{aJe-'^''+'''^+aJeJ^i*-'''^). 

V-configuration 

(45) 

^rE'""^ = Ant = -%a(a<„;ae-^^"* + a.^a+e^^"*) + hg,iaj+e'^''+''^ + aje-'^^^'-^''^). 

Farther we will investigate three-level model concentrated on the internal degrees of freedom 
and neglecting center of mass movement. 

4. The Quantum Dynamics of a Three-Level Atom 

In this section we consider the internal quantum dynamics of a three-level atom in the 
fields of two modes of standing wave resonator. This model is more complicated and less 
studied than the well-known two-level system (see pp, ^ and references cited there). So, 
consider the Schrodinger equation 

(46) ^s^^|^ = i/„,i|vl/(i)> 

for the state vector |^ > in the interaction representation, where the interaction Hamiltonian 
Hint due to H43|) - (|45|l contains the following terms 

(47) H,nt = Hi, = %a(va+e*^»* -t- a^.ae'^"*) -f hg,i&-b+e'^^'+'''^ + aje'^^^*-^'^), 

(48) H,nt = Hit = hga{a;^a+e'''^* + a<„;ae-^^"*) + hgt{a~b+e-''''*+''^ + aje'^^'-''^, 

(49) H^nt = Hl^t = -%a(agia+e-'^^»*-f ag,ae^^"*)-t-%b(<Tre&+e^'^''*+^^^ + ae.6e-*^'*-^^^). 
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We will study the internal dynamics of the three-level system when the interaction Hamil- 
tonian becomes time independent. It is the most convenient way to obtain the time evolution 
operator. 

4.1. Exact resonance in three- level atom. In the case of exact resonance A = the 
interaction Hamiltonians (|47|I - H49|I transform into the following ones 

(50) Hint = -ff^t = hgaia-giO.^ + aigo) + %{,((Tje^^e"''^ + o-^fee^"''^), 

(51) H^nt = Htt = hga{agia+ + fjgia) - hgb{a,M^'^ + aj+e-'^^), 

(52) H,nt = HYnt = -h9a{(yg^a+ + cT;ga+) + hgtiaj+e'^"^ + aje-'^"^). 

Now they do not contain the explicit time dependence and the Schrodinger equation (|46|l 
for H5U|I - H52|) can be integrated formally to give 

where the unitary time evolution operator is defined as 

(53) U{t,to = 0) = exp{-lH^ntt) 

Substituting in H53|l the Hamiltonians (|50|I - H52|I . we obtain the time evolution operators for 
L-, A- and ^-configurations of a three- level atom 

(54) U'^it, to = 0) = cxp[-igat{ag,a+ + &,ga) - ihgi,t{a,,b+e'^'^ + aje-'^'^)], 

(55) U^{t, io = 0) = exp[-iga<('3-gia+ + aigo) ~ igbt{aieb^ e'^"^ + aetbe^^^"^)], 

(56) U^{t, fo - 0) = exp[igat{ag,a+ + ag,a) - ihgi,t{a,,b+ e'^^ + aje-'^"^)]. 
The initial state vector 

(57) 1^(0) = \'^f^eld > (E>\'^atom > = 

cxD oo oo oo / Ce \ 

y^ y^ Wn„WnJCe|e > +Ci\i > +Cg\g >]\n >= X! X! ^"-"^"t. I '^* ) 1"-'^ > I"-'' >' 

ria— Orib— ria—O rit,— 

represents the direct product of the field state vector 



-^9 



oo oo 



(58) 1"^ field >= X! X! "^^^"^ntlna > |"b > 

Ua—O nfa— 

and the atom state vector 

(59) I'^atom >= [Ce\e > +Ci\i > +Cg\g >]. 

It is a convenient form for utilizing the unitary transformations of initial atom-field state 
vector. So, the dynamics of a three-level atom, in the case of exact resonance, is defined by 
the H54() - (|56() time evolution operators. In the next section we calculate these operators and 
utilize the unitary transformations important for applications of quantum computation. 
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5. Calculation of the Time Evolution Operator and Unitary 

Transformations 

5.1. Calculation of the time evolution operator in terms of the exponential ma- 
trix. Based on relations (|5|)- (|10|l let us rewrite the time evolution operators (|54() - (|56|l in the 
matrix form 



(60) 



/ 



U'^^Ht^to = 0) = e 








■-'^"Pgbtb 


gata+ 





goto, 





(61) 



U''^\t,to = 0) = e 








db 



•r'^'Pgbtb+ 




\ 
9ata 
J 



I 

e'^'^qhtb+ 



(62) C/(^)(t,io=0) 

or in the general form as 







where 



U{t,to = 0) = e 
configurations : 



V 



L 







A 

A+ 

B 
B+ 



9ata 

e-'^'figbtb 
e'^fqhtb+ 



-iStp 



Obtb 



I 




-gata 





-gata+ 





B 



A+ 



\ 
A 

; 



gja 
-e^^fgbtb+ 



V 

~gata+ 

-goto, 

e-'^'fgbtb 



-e-'^'figbtb+ e'^fgbtb+ 
corresponding operators from (|60|I - H62|) . including the corresponding scalar factors for con- 



venience. 



5.2. Semiclassical approach. At the given stage we substantially simplify the situation 
by means of simplifying the field operators. Right hand sides of equations H37() include only 
f expiflt type terms, which describe contribution due to free evolution of the field. On the 
other hand it means that the applied fields are sufficiently intense compared with what can 
be contributed by the atomic transitions. The field operators at time t may then be replaced 
by their averages in the initial state [3j. We assume the fields to be in the coherent states 
(the field of laser radiation) |{a„} > and |{/3„} > (|a|,|/3| >> 1), replace the operators a 
{to = 0) and b {to = 0) by a and /3 in (|^ - l|^ and receive the following 
L-configuration 



JSip 





gptp* 



e-'^'^gf,tf3 \ 



?7(^)(i,to = 0) 




gja* 



gata 

; 



A-configuration 



u^^Ht,to^o) 



( ^ 

\ 



-iSif 





gata* 



gptp* \ 
gata 

; 



^-configuration 

f7(^)(t,io = 0) = e 
In general we tiave 



e-'^'fgptp 

e'^'Pgptp* -g^ta* 

—Qnta 



(63) 
whiere 



(64) 



f/(i,<o-0) 
configurations : 



B \ 
B* A 
A* ) 



A 
A* 
B 
B* 



L 

gata 
gata* 

i-'^^gptp 



A 

gata 
gata* 

— iSif 

— i&if 



9f3tP* 

gptp* e^'^^g^tp 



V 

-gata* 
—gata 

-'"'^gfstp 

iSip, 



The exponential matrix (|63|l can be calculated to give 
(65) 

/ AA* + BB'cos(^/BB*+AA 
U{t, io = 0) - 



BB' +AA' 
iB'sin{VBB*+AA') 



iBsin{^BB*+AA* 
VBB*+AA* 



BA{cos{VBB*+AA*)-l) \ 
BB'+AA' 



cos{^/BB* + AA*) 



iAsin{^BB'+AA') 



V 



^/B B'+AA* ^^^yy -L^ ^ I .... J ^BB' +AA* 

A*B*{cos{^yBB''+AA*)-l) iA" sin{^/ BB* +AA*) BB' +AA' cosj^/BB'+AA*) 

BB'+AA^ ^BB*+AA* BB'+AA' 



5.3. Unitary transformations. So, we found the time evolution operator for the three- 
level atom interacting with two modes of field in the coherent states. As the initial state vec- 
tor is represented for the field in the arbitrary state \'^ field >= J2'^ =o S^=o ^naW;„j|na > 
\ni, >, than it is essential to rewrite this state vector for two modes in coherent state. The 
probability that there are n^ (rib) photons at the initial coherent state \a > (|/3 >) is given 
by the Poisson distribution 0], i. e. 



(66) 



< Ua >"" e 



-<na> 



-,W„ 



< Hb >"" £-<""> 



Ua'. Ub'. 

where < n >= |ap (|/3p)- Due to H66|) we can now write the initial state of field (|58|l and 
atom-field system (I57II in more suitable form 

(67) 1*^,;,,, >^ ^ ^ ^^^^^^^1^ ^^^^^^1^ I'^a > \nb >= \a > |/3 >, 

■ria—O nfo— 



I 



Ub 



(68) 



1^(0) >= \'^f^eld > (^l^^atorn > = 



oo oo 



< Ua >"" e' 



< Ub >"■'• e 



"b p-<nb> 



Ua—O nu—O 



Ub 



[ce|e > +Ci\i > +Cg\g >]\na > \nb >- 



OO OC 

EE 

Ua —0 rib— 



< Ha >"" e" 



< Ub >"■'• e 



"i, p-<n.b> 



ribl 



Ci, \aa > \nb >= 



= [ce\e > +Ci\i > +Cg\g >]\a > \P >^ \ c, |a > |/3 > . 

Now applying H65|l to H68|) we obtain the total state vector |^(<) > of a three- level atom in 

the coherent field of two modes for arbitrary time i > to = 

(69) 

(Ce 
c^ I |a > 1/3 > 
Cn 
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5.4. Unitary transformations for L-configuration. For L-configuration H59|) . (|69|) gives 
|*(t) >= f7(^)(i,io = 0)[ce|e > +c,\i > +Cg\g >]\a > |/3 >= U^'^\t,to = 0) I c- ] \a > \f3 > 
and due to (|65|l and H64() we obtain 

uW^t,to^O) = ^--^—^-—x\\ 

, / gl\P? e-'^'^gagbpa 

+cos{tJgl\a\^+gl\(3\^)\ S^a^ + ^g^^ q 

V e'''^gagbP*a* 5^|ap 

^ / e-^'^^gbP 

-sin{tJgl\a\^+glW)\ e-^'^gbP* ga« 

V gaa* 

Tlius, from tiie time evolution operator one can extract three formal parts: time indepen- 
dent, time dependent as cos(G'i) and time dependent as sin(Gi). For the unitary operator 
t/'^^-'(i,to = 0) an initial moment t = we obtain 

/ 1 

(70) UW(t = 0) =010 

V 1 

and 

/ Ce(0)|a>|/3> 

(71) ^(i = 0)= c,(0)|a>|/3> 

V Cg(0)|a>|/3> 

or 

(72) *(t = 0) = ce(0)|e, a, /3 > +c,(0)|i, a, /3 > +Cg\g, a, /3 >, 

where the following notations are introduced |e > |a > |/3 >= |e,a,/3 >, |i > |a > |/5 >= 
\i,a,f3 > and \g > \a > \f3 >= \g,a,f3 > . (|70|l - (|72ll confirm the main feature of the 
interaction representation, that in the absence of interaction the state vectors do not change 

m 

6. Probability Amplitude Method 

In this section, we present a different but equivalent method to solve for the evolution of 
the atom-field system described in interaction picture 

(73) ,n^^^^ ^ H„,t\^{t) > 

and based on the solutions of the probability amplitudes. The atom-field system is described 
by the interaction Hamiltonian H39|l 

(74) H^^^ = %a('Tg»a+e'^"""'^'''^*+<7»gae-'(""-'^''')*-ag,ae-'("»+'^''')*-Ka,ga+e'(^^"+"'»)*-h 

+hgb{a^eb+e'^^^''-'^"'>*+'^'^+ae^be~'^^''-'^"'>*'^'^'^-a^ebe"'^^'+'^"^^^^^ 

where the RWA should be utilized. Further calculations are based on amplitude method of 
QED 0] and are very similar for different configurations. So we perform the calculations 
only for L-configuration and then present the final result for all of them. 
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6.1. L-Conflguration. In the RWA, the Hamiltonian (|74|) . due to configuration (see Fig. 
1), transforms into the f (|47l) - (|49|l ') interaction operators. EspeciaUy for L-configuration 

(75) H,^t = HL = fl9a{(rg^a+e''^^' + a.^ae-'^"*) + ngb{a,,b+ e'^''+'''^ + a^fee-^^"*-'^'^). 
In the exact resonances case (A = ) (|75|) transforms into 

(76) ifj^t = hgaia-giCi^ + aigCi) + hghidieb'^e''^'^ + a^ibe'''^'^). 

Now acting on the wave vector |e,ria,"b > by the operator (|76|l step by step we obtain 

(77) H^„f\e,na,nb >= figbe"""^ ^J Ub + l\e,na,nb >, 



(78) H.^„^\e,na,nb + l>=hga\/n^^^T\g,na + l,ni, + l> +hgbe * '^y/rib + l\e,na,ni, >, 

(79) H^^^f\g,na + l,7if, + 1 >= hgay/ua + l\i,na,nb + 1 >, 

where Ua and Ub represent the number of hfla and ftilf, photons. From H79|l it is also easy 
to see that for ria = and for arbitrary Ub (let us rib + 1 = m) 

(80) i?f„,|g,0,m>=0. 

Now with the aid of obtained amplitudes (|77|I - H8U|) one proceeds to solve the equation of 
motion (|73|l for state vector of total system |^(i) > , i. e., 

oo 

(81) |*(0 >= Yl *e,n„,n,(i)|e,na,n6 > +^^,„„,„,+l(i)|^,na,7^^, + l > + 

n^— 0,nfa— 
+ '^g,na + l,n, + l{t)\9,na + l,nb + 1 > +*g^o,m(i) 1^7 0, m > 

Now substitute l|^ in (|7^ and then projecting the resulting equations onto < e,na,nb\, 
< i,na,nb + l\,< (7,na + l, rib + lj, < g,0,m\ we can obtain the following system of equations: 

(82) ^e,n^,n,{t) = -ifffce"'^'^ V^^;TT*,,„„^„, + i (t) , 

(83) ^i,ri„,nb + l(i) = -igbe''^'^ Vm + l*e,„„,„,(i) - igaVna + l'^g.n^ + l,nt + l{t), 

(84) ^g,n^ + l,n, + l{t) = "iffae"^''^ V^a + l*i,„„ ,„, + ! (t), 

(85) *g,0,m(0 = 0. 

To solve the system of equations H82|l - (|84|l we should define the initial conditions first. 
Under initial conditions the state vector has the following form (see H67f) ') 



oo oc 



|*(t ^ 0) >== |*atom «) \'^f^eld >= [Ce\e > +C.,\i > +Cg\g >] ® J2 Yl ^"o'^nJUa > |nfc >, 

or in notations of H81(l 

oo 

(86) |^(f = 0)>= Y [Wn^.nbCe\e,na,nb> +Wn^,„^ + iCt\i,na,nb + l> + 

71a, "Mb— 

+Wn^ + l^nt + lCg\g,na + l,nh + 1 >] + Wo.mCglff, 0,m > 

where w„„,«6 = Wn^Wn^. 

Comparing H81|) for t ~ and H86|l we obtain the initial conditions 

(87) ^g„^ „^(i = 0) = W„^,„,Ce, 



'^i.n^jit + lit = 0) = W„^,„, + iQ, 
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(89) 



^g,n„ + l,nb + l(i — 0) — ?i'n„ + l,n6 + lCg 



(90) *g.O,™(t = 0) = *3,0,m(t = 0) = U^O.mCg. 

Now, let us solve the system of equations (|5^ - l|M)l under initial conditions defined by ljS7|l - 
()9fl|l . Multiplying l|82(l - (|84|l by e~** and using the Laplace transformation of the probability 
amplitudes the system of equations transforms to the following: 

(91) e-^*^e.n„,„,(i) = -igbVm^TTer'*e-'^'^^,,„^,„,+i{t), 

(92) e-"**,,„„,„,+i(t) = -^igtVrib + le'^'e'^'^ - igaC'^Wna + l*g,„„+i,„,+i(t), 



(93) e ^*'i'e,„„,„,(t) = -igtVna + le ^**.,„„^„,+i(i). 

Now integrating l|M|) - (|^ by t we obtain 



s*^ 



(s) + J56e-""^V'^fc + l*e,n„,„, + l(s) = *e,n.,n,(^ = 0), 



igbe'-^'^Vnb + l*e,7^„,n6 + l(s)+S*e,n„,ni, + l(s)+*56V?^a + l*g,n„ + l,ni, + l ('S) = *-i,na,nb + l (^ = 0), 
igaVna + l*j,„„^„,+l(s) + S*g,„„+i,„,+i(s) = 'I'g,„„+i_„^+i(i = 0). 

Here '^{s) — J^ *(t)e^'**di is the Laplace transform of 'l'(t). This system of equations we 
also can write in the following matrix form: 



U\ ^,,„„,„,+i(s) 



* 



g,na + l,nt + l 



(S) 






where 



s igae-'^f^/^^;^TT 





iga^/ria + 1 S 

Now if we invoke the inverse Laplace transform we obtain the following expression for 
probability amplitudes 

*e,„„,„,(t) 

(94) 



1 /• / '^na-,nb^e 



where 
(7-1 = 



1 



g'+ga("a + l) 



2 



s' + .gaK + i) 



and A = \J g\{na + 1) + ^^("-b + !)■ The poles of (|94|l coincide with roots of the determinant 
\J and equal 0, iA, — iA. From this, due to the Cauchy second integral theorem it follows 



*e,n„,n,(i) 



9l^ 



cos{Xt) + -^ina + 1)(1 - cos(At)) 
A^ 



*»,n„,n. + l(0)- [z^e-'^'^V^I^TTsZ7l(Ai)J *,,„„,„, + ! (0)^ 



^e-^^'^V^I^TTV^irTT(l - cos(At))] *g,„„+i^„,+i(0)} 



^. 



„„,„,+i(i) = { [^e^'^'^V^irrTsm(Ai)] ^e,n„,„,(0) + cos(Ai)] *,,„„,„,+i(0)- 



i^V'^a + Isin(At) 
A 



* 



*. 



,w 



g,n„ + l,n6 + l 



(0)}, 



^e-'^'^VmTTV^^^TTii - cos(At)) 

2 

cos(Ai) + ||(nb + 1)(1 - cos{Xt)) 
A„ 



■ 9a 



i,na,ntiO)~ i-^Vna + lsin{\t) 



* 



i,na.7ib^^ 



(0)- 



* 



g,na+l,ni, + l 



(0) 



Amplitudes for the configurations A and V are obtained in similar manner. 
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The evolution operator 

Ull{t) Ul2{t) Ul3{t) 
U2l{t) U22{t) U23{t) 
Usiit) U32{t) U33(i) 

for the L, A, V configurations has the foUowing form: 
L-configuration, A = ^J gl{na + 1) + 9l{nb + !)• 

u{^{t) = cos{Xt) + ff (n„ + !)(!- cos(Ai)), 
A^ 

uf^{t) = -i^e-''''^VH;~rTsm(At), 
A 

UisW = ^e-'''^V^^;TTV^^^Tl{l - cos{Xt)), 

A 
"22(0 = cos(At), 

A 

31 W = ^e-''^V^^;TTV^^^Tl{l - cosim, 

W32W = -i^e-'^'^^/na + Ism(Ai), 

An 

uf'3(i) = cos(At) + pK + !)(!- cos(Ai)). 



w 



A-configuration, A = \/ g\{na + 1) + 9\^- 



g2 

<i(t) = cos(At) + p(n„ + !)(!- cos(Ai)), 

yA (^) ^ ^^e-*'''^V^s^n(At), 
A 

<3(i) = ^e-^^'^V^V^imCl - cos(At)), 

"2i(*) = «^e*'''^V^sm(At), 
A 

"^2(0 = cos(At), 

'"23(^) = — i — V^a + Ism(At), 
A 

4iW - ^e^''"V^\A^m(l - cos(At)), 
"^2(*) = -«y\/?^a + Ism(Ai), 

2 

u^3(t) = cos(Ai) + ||n6(l - cos(Ai)). 



V-configuration, A = \J g\na + gfinb + 1). 



„2 

M];;(t) = cos(At) + ^K)(l - cos(Ai)), 

uj;2(i) == -i^e-^'''^VH;~n"sm(Ai), 
A 

^^3^ - ^e-*''^V^i;TTV^(l - cos{Xt)), 



u 



V^{t) ^ -i^e'^'^V^^^iPrTsiniXt), 



X 



^22(0 = cos(At), 

'"23(*) = i^\/n^sin{Xt), 

4W = ^e^''^V^^^TTV^(l - cos{Xt)), 



A 
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^32(*) — i-^ \/n^sin{Xt) , 



u^{t) = cos{Xt) + ^{nb + 1)(1 - cos(Ai)). 

Under the assumptions described above we have constructed evolution operators in ex- 
plicit form. Performed analysis has shown that unlike the case of methods of semiclassical 
analysis — where the atom is quantized whereas the field is classical, the completely quantum 
approach shows that in the absence of fields the atom still oscillates and this must be taken 
into account in quantum calculations. Moreover quantum processing may be performed on 
the level of single photons. 
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